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Abstract 

The mutual dependence of spin-dependent conduction and magnetization dynamics of ferromag- 
nets provides the key mechanisms in various spin-dependent phenomena. We compute the response 
of the conduction electron spins to a spatial and time varying magnetization M(r, t) within the 
time-dependent semiclassical transport theory. We show that the induced non-equilibrium conduc- 
tion spin density in turn generates four spin torques acting on the magnetization-with each torque 
playing different roles in magnetization dynamics. By comparing with recent theoretical models, 
we find that one of these torques that has not been previously identified is crucial to consistently 
interpret experimental data on domain wall motion. 

PACS numbers: 75.45. +j, 72.25.Ba, 75.60.Ch 



1 



Recently, there are emerging interests in the interplay between spin-dependent trans- 
port properties and magnetization dynamics of ferromagnets. Giant magnetoresistive effect 
in magnetic multilayers |l| is one of the examples that the spin-transport is controlled by 
magnetization dynamics (or configurations). Spin angular momentum transfer i, or spin 
torque, manifests the magnetization dynamics controlled by spin-polarized conduction elec- 
trons. There are quite a few closely related phenomena reported recently, e.g., enhancement 
of damping parameters due to spin pumping and reaction spin torques [5|, dynamic 

RKKY interaction Q], spin echo and adiabatic spin torques in a domain wall These 
proposed or observed phenomena motivated us to look for a theoretical framework which 
is capable to address the above phenomena on an equal footing. The essence of the above 
phenomena is to recognize two types of electrons: spin- dependent transport is provided by 
electrons at the Fermi level and magnetization dynamics may involve electrons below the 
Fermi sea. While it is impossible to unambiguously separate electrons of transport from 
electrons of magnetization in a real ferromagnet, it has been conventionally modeled via a 
"s-d" Hamiltonian, 

Hsd = —JexS ■ S (1) 

where s and S are the spins of itinerant and localized electrons, and Jex is the exchange 
coupling strength. In this letter, we show that the above simple s-d model in fact captures 
most of the physics on the interplay between spin-polarized transport of itinerant electrons 
and magnetization dynamics of local moments. We will first derive a linear response function 
for the conduction electron spin in the presence of a time and spatially varying local moment, 
and then by using the same s-d model to calculate the spin torque on the magnetization 
dynamics as a result of the induced non-equilibrium conduction electron spin. Among other 
things, we have found four distinct spin torques on the magnetization. Three of them are 
closely related to previously derived torques by using different methods. One of the derived 
torque is new; it describes the mis-tracking between the conduction electron spin and the 
spatially varying local moment. We further show that our formulation can be conveniently 
applied to study magnetization dynamics. An example of domain wall motion is illustrated 
in the end of the paper. 

The dynamics of the conduction electron will be considered separately from that of local 
magnetization. We treat the itinerant spin s as a full quantum mechanical operator whose 
equation of motion is governed by a transport equation, but we approximate S as a classical 



magnetization vector whose dynamics is much slower than that of itinerant spins, i.e., we 
replace S by a classical magnetization M(r, t) 

where |M(r, t)| = Ms is the saturation magnetization. We first determine the induced spin 
density for a given M(r, and then derive the reaction of the induced spin density to the 
magnetization. 

In the present study, the non-equilibrium conduction electrons are generated by applying 
either a DC electric field or a time-dependent magnetic field. While the electric field directly 
generates the charge and spin currents in conducting ferromagnets, the time-dependent 
magnetic field is to drive the magnetization motion that induces a non-equilibrium spin 
density via "s-d" interaction. The conduction electron spin operator satisfies the generalized 
spin continuity equation, 

ds - 1 

g^ + V-J^-[s,Hsd]-rre{s) (3) 

where J is the spin current operator, and rre(s) represents the spin relaxation due to scat- 
tering with impurities, electrons, etc. By defining electron spin density m(r, t) =< s > and 
spin current density J^{r,t) —< J > where <> represents the average over all occupied 
electronic states, e.g., < s >= Tr(ps) where the trace is over all electronic as well as spin 
states, and p is the density operator, one obtains a semiclassical Bloch equation for the 
conduction electron spin density, 

^^^■^=-^^-"'< 

where the commutator in Eq. (3) has been explicitly calculated by utilizing Eq. (2), and we 
have defined Tex — ^/ Jex- 

Next, we separate the induced spin density m into two terms, 

M(r, t) 

m(r, t) = mo(r, t) + 5m(r, t) = '^o^^ \- 5va{v, t) (5) 

where no is the local equilibrium spin density whose direction is parallel to the magnetization. 
The first term in Eq. (5) represents the equilibrium spin density when the conduction electron 
spin relaxes to its equilibrium value at an instantaneous time t. Since the dynamics of the 
magnetization is slow compared to that of conduction electrons, it is reasonable to assume 



the spin of the conduction electrons approximately follows the direction of the local moment, 
known as the adiabatic process. The second term represents the deviation from this adiabatic 
process. Similarly, we write the spin current density as 

Jiv, t) = Jo(r, t) + 5J{v, t) = -(^BP/e)je ® + 5 J(r, t) (6) 

where e is the electron charge, je is the current density, /i^ is the Bohr magneton, and P 
is the spin current polarization of the ferromagnet. Note that the spin current is a tensor 
that consists of two vectors: the charge current and the spin polarization of the current. 
The first term in Eq. (6) is the spin current whose spin polarization is parallel to the local 
magnetization M(r, t). To solve for the non-equilibrium spin density in a closed form, we 
consider the following simplifications. First, we use a simple relaxation time approximation 
to model the relaxation term in Eq. (4), i.e., we write < r(s) >= 5m(r,t)/rsj where r^j is 
the spin-flip relaxation time. The approximation is necessary in order to obtain a simple 
analytic expression. Second, we only consider the linear response of 5m to the electric 
current je and to the time derivative of magnetization dWl/dt. Since 5m is already the first 
order, d5m/dt will be the order of je ■ dlSA/dt or d'^NL/dt'^ and thus it can be discarded. 
Within the semiclassical picture of the transport, the non- adiabatic current density 5J^ is 
related to the non-equilibrium spin density 5m via 5J^ = —DQ'V5m where Dq is the diffusion 
constant. By inserting Eqs. (5) and (6) into (4) and utilizing the above simplification, we 
obtain the closed form for the non-equilibrium spin density 

DoV^5m —5m xM = ^. _ ^ . V M 7 

TexMs Tsf Ms dt eMs 

One immediately realizes that the non-equilibrium spin density is created by two source 
terms on the right side of Eq. (7): one is the time variation and the other is the spatial 
variation of the magnetization. The solution of the above differential equation depends on 
the detail structure of the magnetization vector. Here we assume that the magnetization 
varies slowly in space, i.e., the domain wall width W of the magnetization is much larger than 
the transport length scale defined in the footnote [9] . In this case, the spatial derivation, the 
first term in Eq. (7), can be discarded P]. Then Eq. (7) becomes a simple vector algebraic 
equation and by using the elementary vector manipulation we readily obtain an explicit 
expression for the non-equilibrium spin density 

Tex ( ^nodM. no ^ ^ dM. UbP^ I^bP^^ ^ 

'(8) 



where ^ = Tex/Tsf- The above induced spin density in turn exerts a spin torque on the 
magnetization. From Eq. (2), the torque is T = —{Jex/hMs)yi x m = — ( Jex//iMs)M x 5m. 
By using Eq. (8), we have 

^ +yMx---^Mx[Mx(j,.V)M]-^Mx(j,.V)M 



1 + dt M2 dt eM| ' ; J ^^^2 

(9) 

There are four terms; the first two are from magnetization variation in time and the last 
two in space. Interestingly, the first two terms are independent of the current. The last two 
terms represent the current-driven effect since they are proportional to the current. We now 
discuss the role of each spin torque below. 

The standard Landau-Lifshitz-Gilbert (LLG) equation consists of a precessional term due 
to an effective field and a phenomenological damping term. In addition to these two torques, 
the above torque T is now added to the LLG equation, 

where 7 is the gyromagnetic ratio, Hg// is the effective magnetic field, a is the Gilbert 
damping parameter. We immediately realize that the first term in Eq. (9) is simply to 
renormalize the gyromagnetic ratio while the second term is to renormalize the damping pa- 
rameter. Thus if we introduce an effective gyromagnetic ratio 7' and the damping parameter 
a', 

7' = 7(1 + 7])-^; 7V = 7(a + ^77) 

where we have defined r] = {jiq/ M^) / {1 + ^2), LLG equation remains in the same form. 
We point out that the modification of the gyromagnetic ratio and the damping parameter 
through the present mechanism is rather small in transition metal ferromagnets. For a 
typical ferromagnet (Ni, Co, Fe and their alloys), Jex ~ 1 eV, r^/ ~ lO^^^s, no/M^ ~ 10~^, 
10^^ and thus t] is about 10~^ and ^t] is of the order of 10~^-much smaller than the 
typical damping parameter of the order of 10~^. Therefore, we conclude that the temporal 
spin torque driven by the exchange interaction only slightly modifies the damping parameter 
and can not be identified as a leading mechanism for magnetization damping. 

At this point, we should compare other theories on the spin torque. Tserkovnyak et 
al. proposed an adiabatic spin pumping mechanism to explain the enhancement of 

Gilbert damping parameters. Ho et al suggested a radiation field induced by magnetization 
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precessional motion of magnets j^. Most recently, a similar s-d model in the presence of 
the time-dependent magnetization has been considered The present approach reduces 
to these theories in the simple limit considered for these two terms. In fact, the idea of this 
temporal spin torque had been suggested earlier: when the magnetization varies in time, the 
spin of the conduction electrons tends to follow the direction of the magnetization with a 
time delay given by spin relaxation time; this phenomenon was named as "breathing Fermi 
surface" We are now able to consider this physics of the enhanced damping on the 

equal footing as the current induced spin torques. 

Our main focus here is the spin torque due to the spatially non-uniform magnetization 
vector, the last two terms in Eq. (9). Since the temporal spin torques can be completely 
absorbed by the re-definition of the gyromagnetic ratio and damping constant, we should 
now just ignore them and concentrate on the role of spin torque generated by the non- 
uniform magnetization. We thus write the full equation for the magnetization dynamics 
below 

__ = -7M X Heff + — M X — X M X ^— x ^— (11) 

dt ' Ms dt M2 V dx ) Ms dx ^ ^ 

where we assume the direction of current x-direction (jg = je^x), bj = Pje/iB/eMs(l + ^^), 
and cj = PjefiB^/GMs{l + C,^). Note that bj and cj have the unit of velocity. The term 
has been already proposed by Bazaliy et al. l3| when they consider a ballistic motion of 
conduction electrons in the half-metal materials. Recently Tatara and Kohno also derived 
similar expression jsl. We have seen that this term describes the adiabatic process of the 
non-equilibrium conduction electrons. The "cj" term is completely new; it is related to the 
spatial mis-tracking of spins between conduction electrons and local magnetization. While 
this term is known in the physics of domain wall resistance jl^ . Ilfil llfi| . it also gives rise a 
non- adiabatic spin torque, the last term in Eq. (11). At first sight, one might think that this 
"cj" term may be discarded since it is much smaller than the term {cj/bj = ^ ~ 10^^). 
We will show below that the terminal velocity of a domain wall is independent of the strength 
of , rather it is controlled by this small "cj" term. Thus, experimental analysis on the 
domain wall motion must include this new "cj term. 

To make a concrete prediction on the domain wall dynamics from Eq. (11), we consider a 
Neel wall in a magnetic nanowire whose magnetization vector only depends on the position 
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along the wire, i.e., M = M(a;, t). The effective field entering Eq. (11) is modeled by 

H M 2 A. 

He// = -Jf^^^ + Jp^^^ - 4vrM,e, + Hexte. (12) 

where is the anisotropy field, A is the exchange constant, and AttM^ is the de- 
magnetization field. In the presence of the spin torque, we follow the Walker's prescription 
of the domain wall motion by introducing a trial function M(6', ip) where {6, (p) are polar 



angles in the following form 



(f = (p{t); Intan = ( x — f ?;(r)(ir ) (13) 



2 W{t) V JO 

The first equation assumes that the projection of the magnetization vector in the domain wall 
on the yz plane is independent of the position. The second equation in Eq. (13) postulates 
that the domain wall shape remains a standard Neel-wall form except that the wall width 
W{t) varies with time and the wall moves at velocity v{t). By placing Eqs. (13) and (12) into 
Eq. (11), and by assuming the domain wall width changes slowly as in the Walker's theory, 
we can find two coupled differential equations for determining the domain wall distortion 
parameters (p{t) and W{t). Interestingly, the expression for the velocity of the domain wall 
at the initial application of the current is |l8| 

while the terminal velocity of the domain wall is 

Vt = V{00) = -— r . (15) 

W{oo)a a 

where W{oo) is the terminal wall width that is slightly smaller than the initial Neel wall 
width iy(0). Equations (14) and (15) reveal the different roles played by the adiabatic {bj 
term) and non-adiabatic (cj term) spin torques: the adiabatic torque is most important at 
the initial motion of the wall while the non-adiabatic cj controls the terminal velocity of 
the domain wall. The adiabatic torque causes the domain wall distortion. The distorted 
domain wall is able to completely absorb the adiabatic spin angular momentum so that the 
net effect of the adiabatic torque on the domain wall velocity becomes null, i.e., domain 
wall stops. In contrast, the non-adiabatic spin torque behaves as a non- uniform magnetic 
field cjd'M./dx that can sustain a steady state wall motion. Although the magnitude of the 
non-adiabatic torque cj is about two orders of magnitude smaller than adiabatic torque 6j, 



the terminal velocity is inversely proportional to the damping parameter which makes the 
velocity comparable to bj. 

Finally, we emphasize that the present study has resolved an outstanding mystery between 
the recent experimental observation and the theoretical prediction based on the adiabatic 
spin torque. It has been recognized that a critical current density of the order of 10^ — 



10^^ A/ cm'^ is required to move a perfect domain wall if we only use the adiabatic spin 

torque bj. Experimentally, a velocity about 3 m/s was observed in a NiFe nanowire when a 
current density 1.2 x lO^A/cm^ was applied. This velocity had been assumed to relate with 
bj 13] in spite of the apparent qualitative and quantitative disagreement between theory and 
experiment. Here, we have pointed out that bj is simply an initial velocity of the domain 
wall and the measured velocity was the terminal velocity. For the experimental current 
density of 1.2 x lO^A/cm'^, the adiabatic spin torque alone is unable to sustain a constant 
velocity. By including a small non- adiabatic torque cj, we find the domain wall velocity is 
now cj/a in the absence of the magnetic field, see Eq. (15). Although the numerical values 
of both the exchange constant Jex and the damping parameter a are not precisely known in 
ferromagnets, we estimate that the wall velocity should be 6 ~ 60 (m/s) for the above current 
density if we use the parameters indicated before (taking a = 0.01 ~ 0.1 for permalloy). 
While the experimental velocity is smaller than our estimated value, it is reasonable that 
we do not include any defects that may reduce the observed velocity significantly. 
The research was supported by NSF grants ECS-0223568 and DMR-0314456. 
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